Methods for Parallel Integration of Stiff Systems of ODEs
نویسنده
چکیده
This paper presents a class of parallel numerical integration methods for stiff systems of ordinary differential equations which can be partitioned into loosely coupled sub-systems. The formulas are called decoupled backward differentiation formulas, and they are derived from the classical formulas by restricting the implicit part to the diagonal sub-system. With one or several sub-systems allocated to each processor, information only has to be exchanged after completion of a step but not during the solution of the nonlinear algebraic equations. The main emphasis is on the formula of order 1, the decoupled implicit Euler formula. It is proved that this formula even for a wide range of multirate formulations has an asymptotic global error expansion permitting extrapolation. Besides, sufficient conditions for absolute stability are presented.
منابع مشابه
On second derivative 3-stage Hermite--Birkhoff--Obrechkoff methods for stiff ODEs: A-stable up to order 10 with variable stepsize
Variable-step (VS) second derivative $k$-step $3$-stage Hermite--Birkhoff--Obrechkoff (HBO) methods of order $p=(k+3)$, denoted by HBO$(p)$ are constructed as a combination of linear $k$-step methods of order $(p-2)$ and a second derivative two-step diagonally implicit $3$-stage Hermite--Birkhoff method of order 5 (DIHB5) for solving stiff ordinary differential equations. The main reason for co...
متن کاملImplementation of four-point fully implicit block method for solving ordinary differential equations
This paper describes the development of a four-point fully implicit block method for solving first order ordinary differential equations (ODEs) using variable step size. This method will estimate the solutions of initial value problems (IVPs) at four points simultaneously. The method developed is suitable for the numerical integration of non-stiff and mildly stiff differential systems. The perf...
متن کاملAdaptive Linear Equation Solvers in Codes for Large Stiff Systems of ODEs
Iterative linear equation solvers have been shown to be effective in codes for large systems of stiff initial-value problems for ordinary differential equations (ODEs). While preconditioned iterative methods are required in general for efficiency and robustness, unpreconditioned methods may be cheaper over some ranges of the interval of integration. In this paper, we develop a strategy for swit...
متن کاملLinearly Implicit Discrete Event Methods for Stiff ODEs. Part I: Theory
This paper introduces two new numerical methods for integration of stiff ordinary differential equations. Following the idea of quantization based integration, i.e., replacing the time discretization by state quantization, the new methods perform first and second order backward approximations allowing to simulate stiff systems. It is shown that the new algorithms satisfy the same theoretical pr...
متن کاملAn Almost L-Stable BDF-type Method for the Numerical Solution of Stiff ODEs Arising from the Method of Lines
A new BDF-type scheme is proposed for the numerical integration of the system of ordinary differential equations that arises in the Method of Lines solution of time-dependent partial differential equations. This system is usually stiff, so it is desirable for the numerical method to solve it to have good properties concerning stability. The method proposed in this article is almost L-stable and...
متن کامل